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ABSTRACT 

We use the ultra-deep HUDF09 and the deep ERS data from the HST WFC3/IR camera, along 
with the wide area Spitzer /IRAC data from GOODS-S to derive SEDs of star- forming galaxies from 
the rest-frame UV to the optical over a wide luminosity range (M1500 ~ —21 to M1500 ~ —18) from 
z ~ 7 to z ~ 4. The sample contains ~ 400 z ^ ^ 120 z ^ 5, ^ 60 z ^ 6, and 36 prior 
z ^ 7 galaxies. Median stacking enables the first comprehensive study of very faint high-z galaxies at 
multiple redshifts (e.g., [3.6] = 27.4 ±0.1 AB mag for the M1500 ~ -18 sources at z ~ 4). At z ^ 4 
our faint median-stacked SEDs reach to ~ 0.06L*^4 and are combined with recently published results 
at high luminosity L > L* that extend to M1500 ~ —23. We use the observed SEDs and template 
fits to derive rest frame UV- to-optical colors {U — V) at all redshifts and luminosities. We find that 
this color does not vary significantly with redshift at a fixed luminosity. The UV-to-optical color does 
show a weak trend with luminosity, becoming redder at higher luminosities. This is most likely due 
to dust. At z > 5 we find blue colors [3.6] — [4.5] ~ —0.3 mag that are most likely due to rest-frame 
optical emission lines contributing to the flux in the IRAC filter bandpasses. The scatter across our 
derived SEDs remains substantial, but the results are most consistent with a lack of any evolution in 
the SEDs with redshift at a given luminosity. The similarity of the SEDs suggests a self-similar mode 
of evolution over a timespan from 0.7 Gyr to 1.5 Gyr that encompasses very substantial growth in the 
stellar mass density in the universe (from ~ 4 x 10^ to ~ 2 x 10^ MqMpc"^) 

Subject headings: galaxies: evolution — galaxies: high-redshift 



1. INTRODUCTION 

As a result of the installation of the WFC3/IR cam- 
era on the Hubble Space Telescope, it is now possi- 
ble to obtain extremely deep, high resolution observa- 
tions in the near infrared. These observations match 
the exquisite optical data provided by the ACS cam- 
era over deep and ult ra- deep field s like the GOODS an d 
the HUDF (Giavalisco et al.ll2004l : fBouwens et al.ll2007[ ). 
The upgrade has enabled several teams to extend their 
searches for Lyman Break Galaxies (LBGs) to very high 
redshifts (e.g., Bouwens et al. 2010b; Oesch et al. 2010, 
[2OII; McLure et al. 2010; Yan et al. 2010|; |Bunker et alj 
[2OIO; Finkelstein et al. 2010). In combination, these two 
cameras have allowed for high S/N measurements of the 
rest-frame UV light of high redshift galaxies over a wide 
wavelength baseline. 

As a result of these new observations, the UV LFs 
at 2: > 4 are known to high accuracy. Another direct 
result from these observations is the determination of 
the UV slope of the SEDs, genera lly characterized by 
the parameter ^ (f x ex A^, e.g., iMeurer et al.l 119991 : 
iBouwens et al.|[2QQ9l ). This slope has been shown to be- 
come steeper (bluer colors) both with decreasing lumi- 
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on the dust reddening so this has been interpreted as an 
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increase in the dust content of galaxies w ith cosmic time 
and as galaxies increase in luminosity (jBouwens et al.l 
IMS [Finkelstein et al. 2011). 

Since the rest-frame UV light primarily probes young 
short-lived stellar populations (mostly O and B stars), 
it is interesting to extend these observations to longer 
wavelengths and study the relative importance of more 
evolved stellar populations. This is very challenging 
though and is now only possible thanks to the incredi- 
ble sensitivity of the IRAC camera on the Spitzer Space 
Telescope. Some of the first observations of rest-frame 
optical light from z > A galaxie s showed qui t e substan- 
tial U V-to-optical colors (Evle s et al.l 120051 : lYan et aTl 
'2005). Under the usual assumption of smooth, ex- 
ponentially declining star formation histories (SFHs), 
these colors were interpreted as indicative of an evolved 
stellar population (ages > 100 Myr) in combination 
with the younger populations probed by the UV-light. 
Follow-up work on z ^ 4 — 7 galaxies continued to re- 
veal systems wi th nioderate size UV-to-optical breaks 
(e.g.lEy les et ar'200?;"Yan et al."2006VLabbe et al.""2006"; 
iStark et al. 2009: G onzalez et al. 201 0), even possibly to 
z ^ (|Labbe et al . 2010a). However interesting, these 
studies are only possible for the brightest z > A galaxies, 
due to the extreme depths required. 

Some alternative explanations have been proposed for 
the red UV-to-optical colors observed in z > 6 galaxies. 
For example, for young star-froming galaxies, one might 
expect strong line emission. At the right redshift, these 
lines could contaminate the 3.6 and 4.5 jim IRAC fluxes, 
resulting in red UV-to-optical colors similar to those pro- 
duced by the Balmer break (Schaerer & de Barros 20091 
l2010f ). For emission lines to be strong enough to domi- 
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nate the rest-frame optical fluxes, extremely young ages 
(~ 6Myr) are required. Such models have ages shorter 
than the dynamical times of these high— z sources and 
are also somewhat uncertain. It would be surprising if 
the majority of these galaxies, which are selected over 
a redshift range that spans > 200 Myr, had such young 
ages. Even if rest-frame optical emission nebular emis- 
sion lines do not dominate the measured IRAC fluxes, 
it is expected that they will contribute some of the flux 
measured in the mid-IR, but the exact contribution cur- 
rently remains hard to establish at these high redshifts. 

It has been routine, then, to compare the colors of these 
high redshift galaxies with synthetic stellar population 
(SSP) models of relatively evolved galaxies with smooth 
SFHs to derive physical properties from the observations. 
By studying LBG samples at different redshifts in the 
aforementioned way, one intriguing result has emerged. 
At a given intrinsic luminosity, the specific Star Forma- 
tion Rate (SFR) of gal axies seems to remain remarkably 
constant w i th redshift ([Stark et aT]l2009l : iGonzalez et al.l 
120101 . 120111 : iMcLure et all 120111 : slightlv modified to ac- 
count for the latest dust corrections by "Bouwens et alj 
[2011). This result has called the attention of theorists 
and has proved very hard to reproduce in simul ations 
fe.g.. lWeinmann et al.llMllKhochfar fc Silkll2011f ). The 
assumptions used to derive these result vary somewhat 
from group to group. In this work, we present the ba- 
sic observations that lead to such a result, namely the 
observed SEDs of galaxies at z ~ 4 — 7. 

In the present work we study the observational prop- 
erties of the typical high redshift galaxy. We do this 
by splitting our large sample in both redshift and UV- 
luminosity bins and constructing mean SEDs. The large 
number of galaxies allows us to study the typical col- 
ors at high S/N but also extend our measurements to 
faint limits that are currently inaccessible on an individ- 
ual basis with the current data. Although the scatter 
about the typical properties is also highly interesting, we 
defer their study to a future work which relies on even 
deeper Spitzer observations. In section 2 we present the 
Data and in section 3 the sample selection and photom- 
etry. Section 4 describes the stacking procedure and we 
discuss the stacked SEDs and colors in section 5. A sum- 
mary is presented in section 6. 

Throughout, we use a (i^o, ^m, ^a) = (70km s~^, 
0.3, 0.7) cosmology when necessary and we quote all mag- 
nitudes in the AB system (Oke & Gunn.1983.) . 

2. DATA 

In this paper, we make use of the large number of 
z > 4 galaxies found in the ultra-deep HUDF data and 
wide-area ERS data to construct median SEDs for star- 
forming galaxies at z ~ 4 — 7. The deep optical, near-IR, 
and deep IRAC coverage over these fields allow us to 
create SEDs that cover both the rest-frame UV and the 
rest-frame optical down to very faint fluxes. 

2.1. HST ACS and WFC3/IR 

In constructing our median SEDs, we make use of 
the deep optical and near-IR HST data from ACS and 
WFC3/IR. Over the ERS field, both the ACS optical 
(^435^^606^775 ^850 ) and the WFC3/IR (Igg Ji25^i6o) im- 
ages from HST reach depths of ~ 28 mag (5 <j measured 
on 0.35''-diameter apertures). In the HUDF field, the 



TABLE 1 
Sample Summary. 





z - 4 


z ~ 5 


z ~ 6 


2-7^ 


ERS 


270 / 524 


77 / 123 


20 / 32 


15 / 18 


HUDF 


137 / 205 


41 / 55 


38 / 51 


21 


TOTAL 


407 / 729 


118 / 178 


58 / 83 


36 



Note. — Summary of the number of sources in the 
sample. IRAC photometry of these sources requires 
fitting and subtraction of the flux from surrounding 
foreground neighbors. This is not possible in all the 
cases. The table shows the number of sources with 
clean IRAC photometry / total number of sources 

found in the fleld. 

^ We adopt the sample presented in ILabbe et al.l 
[2010^. 

ACS optical data (BViz) are 1.5 — 2 mags deeper, and 
the WFC/IR data (Yios -^125^160) are 1.5 times deeper 
than the ERS. The HUDF was observed with the Yios 
filter instead of the Ygs filter used in the ERS. 

2.2. Spitzer /IRAC 

Because these fields are located inside the GOODS-S 
field, they both have deep Spitzer /lYl AC coverage that 
amounts to ~ 23.3 hr over the ERS and twice that for the 
HUDF which was observed to the same depths but in two 
different epochs with the IRAC camera rotated by 180°. 
For IRAC channel 1 images (Ao = 3.6/im), each single 
epoch reaches depths of 27.8 mags, measured over 2.5'^- 
diameter apertures (Icr). For channel 2 (Ao = 4.5/im) 
this limit corresponds to 27.2 mags. 

3. SAMPLE SELECTION AND BASIC PHOTOMETRY 

The criteria used to identify high redshift star form- 
ing galaxies has already been presented in numerous 
previous works. For details on the selection procedure 
and discussion of sample contamination please see 
I Bouwe ns et al.' (2007*, 2010a). In short, the selection cri- 
teria used to find the sources presented here is as follows: 

z A 5-dropouts: 

(^435 - ^606 > 1.1) A [^435 - 1^606 > (^606 - ^85o) + 1-1] 
A (1^606 - ^850 < 1.6) 

z ~ 5 F-dropouts: 

{[V^606 - ^775 > 0.9(^775 - ^85o)] V (^606 - ^775 > 2)} 
A (1^606 - ^775 > 1.2) A (^775 - ^850 < 1.3) 

z ~ 6 z-dropouts: 

(^775 - ^850 > 1.3) A (Z850 - -^125 < 0.8) 

The combined samples contain a total of 729 sources 
at z ~ 4, 178 at 2: ~ 5, and 83 at z ~ 6 (see table [1]). 

We complement the 2: ~ 4, 5, and 6 sam ples with the 
sample of 36 z ~ 7 sources presented by Labbe e t al.l 
(2010a). We have adopted the z ~ 7 stacks presented in 
that work. 

At these depths, the Spitzer /lH AC images are fairly 
crowded due to the size and shape of its PSF. As a con- 
sequence, standard aperture photometry of the sources 
generally results in fluxes that are contaminated by flux 
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Fig. 1. — The observed magnitude distribution of the sample. The reference filters for each redshift sub-sample have been chosen to be 
close to rest-frame A = 1500 A. For the z ^ 6 sources in the ERS we used the Yqs filter whereas the Yios filter was used for the sources in the 
HUDF. The open histograms include all the sources and the gray filled histograms include only the sources with clean IRAC photometry. 
This selection does not introduce biases relevant to the determination of the mean rest-frame optical colors of these sources. 



coming from umelated sources that are nearby in the 
sky. We are able to de-blend the fluxes from the differ- 
ent sources by making use of the higher resolution infor- 
mation available through the HST images. Our method 
has been explained in previous works by our group (e.g 
iLabbe et al.ll2QQll2QlQbl : IGonzalez etlIl[Mll2QTl[ ) but 
we give a short description here. We start by regis- 
tering the IRAC images to the higher resolution HST 
image that will be used as a light profile template. In 
this case we use the WFC3/IR i^ieo-band image as our 
template because it is the closest in wavelength and has 
high signal to noise. W e use the SExtractor software 
(jBertin fc ArnoutsI Il996f ) to create segmentation maps. 
We isolate all the sources in a small area around each 
galaxy in our sample using these segmentation maps. We 
then convolve them with a kernel derived from the PSFs 
of the HST and the IRAC images. We fit for the normal- 
ization fluxes of all sources (including the source we are 
trying to clean) simultaneously. Finally we remove the 
flux from the unrelated sources to obtain a clean image 
of the galaxies in our sample. Later we use these cleaned 
images to produce median stacks. 

The procedure described above does not always pro- 
duce acceptable results. The reasons why it can fail 
vary but the most common is the presence of a source 
that is too bright, extended, and close to our region of 
interest. In these cases the amount of flux subtracted 
off our source of interest is large and more uncertain 
and in most cases results in large residuals. We have 
carefully inspected the residual images (all sources sub- 
tracted, including the one in our sample) and the cleaned 



images (only unrelated neighbors subtracted) of each of 
the sources in the original sample and come up with a 
sub-sample of sources with clean IRAC photometry. This 
"by eye" inspection criterion generally results in good 
agreement with a criterion that selects only the images 
with the lowest residuals but allows us to also re- 
move sources that, even though have moderate val- 
ues, present mild to strong residuals located very close 
to the source. Overall, this reduces the sample to 60% 
of the original (see Table [1]). Since the criterion adopted 
here depends on the complexity of the unrelated sources, 
it is not expected that it should introduce any biases 
that would be relevant in the determination of the mean 
rest-frame optical colors (see Figure [T]). 

4. STACKED PHOTOMETRY AND MEDIAN SEDS 

It has been shown that in star forming galaxies 
the total SFR appears to correlate with stellar mass, 
the so called main sequence of star forming galaxies 
(NoeskieFaDH^. This sequence has been observed 
at high-redshifts {z > 4) but contrary to what has been 
observed at z < 2, the normalization of the relation at 
fixed stellar mass does not seem to evolve over 2 < z < 7 
(|Stark et al.l I2QQ9I : iGonzalez etld] I2Q11I) . The lack of 
deeper rest-frame optical from Spit zer /IRAC for these 
sources is one of the big Hmitations to study this sequence 
to fainter magnitudes. In terms of the bulk of the popu- 
lation, however, progress can be made by stacking large 
numbers of sources to study their mean properties. In 
this section we seek to construct the median SEDs of 
galaxies at z ^ 4. 

Due to the nature of the dropout search and the filters 
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available, our full sample is naturally divided in redshift 
bins cen tered at {z) = 3.8, 5.0 , 5.9 plus the {z) = 6.9 sam- 
ple from lLabbe et al.l (|2QlQa[ ). In this section, we further 
split each redshift sub-sample based on their observed 
rest-frame UV- luminosity. The filters used as reference 
UV- luminosity are the ^775 filter for the z ^ A sources; 
the ^850 for the z ~ 5 sources; the Igg for the z ~ 6 
sources in the ERS and the Yios for the z ~ 6 sources in 
the HUDF; and finally, the J125 filter the z ~ 7 sources. 
These filters were chosen to be close to rest-frame 1500 A 
at the different redshifts and make the splitting criterion 
fairly uniform across redshifts. These bins are presented 
in Figure [H 

In the case of the HST ACS and WFC3 imaging, the 
fluxes utilized correspond to the MAG_AUTO values as 
determined using the SExtr actor code. For each UV- 
luminosity bin, the median flux in each band was de- 
termined from the individual flux measurements of the 
sources in the bin that also have clean IRAC determi- 
nations. The uncertainty in the median was estimated 
through a bootstrap process. In each realization, Ni^in 
fluxes are drawn with replacement from the Ni)in values 
available and perturbed by their individual uncertain- 
ties. Generally the individual HST flux determinations 
have very high signal to noise. There are however, sys- 
tematic uncertainties in the absolute calibration of the 
different bands and instruments across the SED. These 
uncertainties amount to at least 5% which we adopt as 
our minimum allowed uncertainty in any given band. 

In the case of the IRAC imaging, most sources are 
too faint to be individually detected even in these deep 
^ 23.3 hr integration images. Instead of using individ- 
ual fluxes, we first median combine (pixel by pixel) the 
cleaned stamps of the sources within a bin that have 
clean residuals from our IRAC cleaning procedure. We 
perform standard aperture photometry on the stacked 
images using a circular aperture of 2.5''-diameter and 
correct the fluxes to total assuming PSF profiles. This 
amounts to a factor 1.8 in both the [3.6] and the [4.5] 
IRAC filters. The local background is estimated from 
the stacks, in a ring between 6''and lO^'diameters around 
the stamp center, where the source is expected to be, and 
then subtracted off the image. 

To estimate the uncertainty in the IRAC stacks, we 
also create bootstrap realizations. We treat the two 
epochs of IRAC data over the HUDF as independent 
stamps in the stack, thereby effectively weighting with 
exposure time, i.e., Nun = Ners ^ '^Nhudf- In 
each bootstrap draw, then, we draw with replacement 
Ners + "^Nhudf such stamps, and estimate fluxes as 
described above. The image noise in each stacks is ex- 
pected to decrease with the total equivalent exposure 
time approximately as 1/ -sjt^^x'p.eq {texp,eq — {Ners + 
"^Nhudf) X 23.3 hr). The actual photometric uncertain- 
ties that we derive are always greater than this because 
they also include the intrinsic scatter in the fluxes of the 
sources been stacked together (see Appendix). Stamps 
for all the IRAC median stacks along with the derived 
magnitudes and uncertainties are presented in Figures [21 
[21 and [H Table [2] is a c ompilation of all the stacked pho- 
tometry including the iLabbe et al.l (|2010a[ ) stacks and 
information on the number of sources that go in each 
stack. It can be seen in the Figures that even for sources 



as faint as 27 mag, the stacks show good detections in 
both IRAC bands. 
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[3.6] 



[4.5] 



27 ERS + UDF 
'~629.1 eq. hrs 



^775 ^26 

87 ERS + 12 UDF 
~2586.3 eq. hrs 



^775 ^27 

107 ERS + 25 UDF 
~3658.1 eq. hrs 



^775 ""28 

45 ERS + 43 UDF 
~3052.3 eq. hrs 




24.20 -±0.11 



|.^24.23±0 .17j 



25.37 ±0.07 



25.32 ±0.07,1 



♦ 

26.63 ±0.( 



26.87±0.1 



27.40 ±0.13 



27.34 ±0.18 



Fig. 2. — Median stacked images of the 2; ~ 4 sources. The sam- 
ple has been spht according to their irrs^AB magnitude (approx. 
rest frame 1500 A) in bins of 1 mag. The stamps are lO^^on a side. 
The number of sources stacked along with the total equivalent in- 
tegration time has been included to the left of the stamps. The 
magnitudes included in the stamps were measured using circular 
apertures of 2. 5^^ in diameter and corrected to total assuming stel- 
lar profiles which amounts to ~ —0.6 mags in both the [3.6] and 
the [4.5] channels. The background was estimated from the stacks 
and subtracted off the images. To estimate the uncertainties on the 
fluxes measured we created 200 random realizations of the stack. 
In each realization, sources were drawn with replacement and then 
median stacked pixel by pixel. 



5. REST-FRAME UV-TO-OPTICAL COLORS OF HIGH-Z 
STAR-FORMING GALAXIES. 

The observed colors of galaxies can give us important 
information about the stellar populations that compose 
them. Usually, this is accomplished by comparing their 
colors with those of Synthetic Stellar Population (SSP) 
models. In this work, we chiefly explore the observed UV- 
to-optical colors and what can be learned from them and 
from their dependence on luminosity and redshift. When 
models are invoked, we use the Bruzual & Chariot! (I2QQ3L 
BC03) models and we assume a Salpeter (1955) IIVIF, 
constant SFH, and 0.2 Zq metallicity, as these quantities 
cannot be constrained from our dat a. To derive best fits 
we use the SED fitting code FAST fiKriek et al.lf2QQ9> ). 

The rest-frame UV-to-optical color is of particular 
interest at high redshifts. Several studies have de- 



0.5 ma^) in 2; > 4 galaxies (IEv 


ies et all I2005L 


2007; 


Yan et al.ll2005UStark et al.l "20091: 
20111: ILabbe et al. 2006. 2010a: McI 


Gonzalez et al.l 
jure et al.ll2011h 


2010, 
. The 



usual interpretation is that these colors derive from large 
Balmer breaks which are associated with evolved stel- 
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24.36 ±0.10 



[4.5] 

24.70 ±0.13 



^26 



^850 

18 ERS + 2 UDF 
~512.6 eq. hrs 



25.12 ±0.13 
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^850 ^ 



27 



43 ERS + 8 UDF 
~1374.7 eq. hrs 



# 

:a§.28±0.11 



26.56 ±0.1^ 



Zo 



^28 



^850 

11 ERS + 9 UDF 
~675.7 eg. /irs 



27.03 ±0.29 



27.62 ±0.50 



Fig. 3. — As in Figure[2]but for the z ~ 5 sample. In this case the 
sample was split according to their zsso magnitude to approximate 
rest-frame 1500 A. 
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3 ERS + 3 UDF 
~209.7 eq. hrs 
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11 ERS + 2 f/DF 
~349.5 eq. hrs 



y-28 

4 + 6 UDF 
~372.8 eg. hrs 



[3.6] 

25.3a jt0,22 



[4.5] 



25.76 ±0.24 




26.32 ±0.31 



,27.44 ±0.54 



Fig. 4. — As in Figure [2] but for the z ~ 6 sample. In this case 
the sample was split according to their Ygg magnitude in the case 
of the ERS sources and lios for the HUDF sources (approximately 
rest-frame 1500 A). 

lar populations with ages > 200 Myr. This is surpris- 
ing given how young the universe is at the time these 
galaxies are being observed. The main limitation for 
field studies is the depth of the Spitzer /IUAC imaging 
available, which constrains these studies to the brightest 
LBGs. Searches behind clusters provide an alter native to 
probe the fainter populations (see, e.^., lEllis~et aL 200 li 
Bradlev et all 12008: "Richard et alj 120081 : iLaporte et aLl 
Ml and the recent tentative determination of a 1.5 mag 
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Fig. 5. — Stacked SEDs of z ~ 4 galaxies in 7 bins of UV- 
luminosity. The faintest SEDs (blue solid symbols and lines) cor- 
respond to the SEDs derived from our sample. This sample has 
a mean redshift of {z) = 3.8 ± 0.3 and was split in bins of ob- 
served Az775 = 1 mag, approximately covering the UV-luminosity 
range —21 < M1500 % —18. The brightest stacks (gray open sym- 
bols and dashed lines) correspond to the median SEDs presented 
in [Lee et al . (2011), for a sample of galaxies with mean redshift 
{z) = 3.7 ± 0.4. These SEDs allow us to extend the l uminosity 
range to brighter luminosities up to M1500 ~ —23. In lLee et al.l 
(|201in the sources were split in six bins according to their observed 
/-band luminosity. For clarity, we only show every other bin but 
all are considered in Figure [21 The x-axis shows wavelength and 
the approximate reference filter in their study. Best fit BC03 mod- 
els are included for references (the H and Xg- band fiuxes were 
excluded from the fits because they are likely biased - see text, sec- 
tion 5.1). A trend of bluer UV-slopes and smaller Balmer breaks 
at fainter luminosities is observed. 



Balmer break in a z 6 galaxy by Richard et al.ll20l"TI ) 
but the number of sources amenable to such studies, re- 
main low. We work around this limitation by stacking 
large numbers of sources which allows us to estimate the 
typical rest-frame UV-to-optical colors of faint z > 4 
sources. 

We start by presenting the results fou nd at z ^ 4 where 
the similar and complementary work of iLee et al.l (|2011f ) 
allows us to study this colors over an unprecedentedly 
large range of UV-luminosities. 

5.1. The UV-to-Optical Colors at z ^ 4. 

Our sample of 5— dropouts has a mean redshift of 
(z) = 3.8 ±0.3. We divide this sample based on their ^775 
luminosity (which approximately corresponds to rest- 
frame 1500 A) in four bins of Ai^^^ = 1 mag. These 
are the same bins as shown in the histogram in Figure [TJ 
The bins at ^775 > 28 mag are too faint to be detected in 
Spitzer /IRAC and will not be included in the analysis. 
Assuming the mean redshift of the sample, these bins cor- 
respond approximately to M1500 = —21,— 20, —19, and 
— 18, which covers the moderately brigs: ht to faint end 
of the population (M^^qq = —20.24, iBouwens et al.l 

I2QQ3). , [ ^ 

In a recent study, ILee et al.l (|2011l ) focus on the me- 



6 



dian SEDs of the brightest z ^ 4 galaxies found in the 
ground based NOAO Deep Wide-Field Survey. Their 
sample has a mean redshift (z) = 3.7 ± 0.4, and their 
SEDs also sample the rest-frame UV and optical using a 
combination of optical filters: Bw^ I (Mosaic Cam- 
era), near-IR filters: J, H, Ks (NEWFIRM Camera), 
and the Spitzer/IRAC mid-IR channels 1 through 43- 
Similar to the present work, they divide their sample ac- 
cording to the /-band luminosity, although with uneven 
bins. Their sources are intrinsically brighter than in this 
work, corresponding to L > I/*=4 sources roughly cover- 
ing the range -23 < Mi,n n < -21. 

The combination of the TeeetaD (l 2Qllf ) sample with 
our data spans the UV-luminosity range: —23 < M1500 ^ 
— 17.5. This represents an unprecedented dataset to 
study the SEDs, and especially the UV-to-optical col- 
ors, of z ~ 4 galaxies over a large range of luminosities. 
The SEDs from both works are presented in Figure [H In 
this figure the i/ieo is shown with open symbols since it 
is likely biased and should be ignored as is explained in 
detail later. An overall trend to bluer colors (both in the 
rest-frame UV and the rest-frame UV-to-optical colors) 
towards fainter luminosities is already apparent in this 
figure. 

We quantify the UV- to- Optical colors by estimating 
the interpolated rest-frame {Un — Vn) colors for the 
median-stacked SEDs in both sets. Figure [6] shows a not 
very steep but systematic trend of redder {Un — Vn)rest 
colors as a function of increasing UV-luminosity. 

The interpolated {Un — Vn) colors were determined 
from the best fits to the full observed SED, where the fits 
correspond to a linear combination of a set of template 
SEDs. These SEDs correspond to the default template 
set from the photometric code EAzY (Brammer et al. 
I2OO8). This template set was constructed and calibrated 
empirically to determine accurate photometric redshifts 
of galaxies of various types, from very blue starbursts to 
red ellipti cals and has been s hown to work well up to 
z - 5 - 6 f Brammer et al.l[2QQ8 ). For the Un and Vn fil- 
ters we have assumed narrow ideal filters that correspond 
to step functions of width 100 A, centered at 3500 A and 
5500 A for the Un and Vn filters respectively. The colors 
were calculated from the best fit combination of tem- 
plates using the filters described above. To estimate the 
uncertainty in the {Un — Vn)rest estimates, we perturbed 
the photometry and the redshift of each median-stacked 
SED within the allowed uncertainty, and a new fit was 
obtained. A new set of colors was measured from this 
best fit and the procedure was repeated a large number 
of times to derive the confidence intervals. A similar pro- 
cedure was used to derive rest-frame M1500 magnitudes 
for an ideal filter centered at 1500 A and with a width of 
100 A. 

For the data-point at M1500 ~ —19, only an upper limit 
can be obtained due to the limitations of the template set 
used here. The EAzY default template set only includes 

^ Although relatively deep IRAC channels 3 &; 4 data (sampling 
5.8 jim and 8.0 jim respectively) is available for our sample in the 
GOODS-S, this data is still > 0.5 mag shallower than the channel 
1 &; 2. It is, therefore, much harder to obtain reliable stacked 
photometry in those bands, especially considering that our sample 
is much smaller and intrinsically fainter than the Lee et al. (2011) 
sample. We have not included these bands in the median-stacked 
SEDs presented in this work. 



SEDs with {Un - Vn)rest > 0.3. The SED at M1500 - 
— 19 is likely to have a {Un — Vn)rest color of ~ 0.3 mag, 
so no reliable lower limit can be obtained. We choose 
to show a 2 cr upper limit. Other more sophisticated 
algorithms to determine interpolated rest-frame colors 
for an observed SED y ield very similar results (e.g., the 
InterRest code jTavlor e t al. 200 ^ which imp lements the 
algorithm described in Rudnic k et al. 1 120031 ). However, 
such algorithms also depend on the template set resulting 
in similar limitations. 

It is important to note that, when deriving this color 
from best fits, we have ignored the ij' iBn—band in our 
stacks and the H and i^Tg— bands in the lLee et al.l (|2011f ) 
stacks. Given the redshift distribution of the sample it 
is expected that for some fraction of the sources (the 
ones at the lowest redshifts) the i^— band fluxes mea- 
sured will include some rest-frame optical light. As a 
result, the estimated median rest-frame UV flux, will be 
biased towards larger values. To attempt a correction for 
this would require us knowing the exact redshift distribu- 
tion of the sample, but even then a correction would be 
highly uncertain, so we opt to remove these photometric 
points from all fits^ 

In recent works jBouwens et al.l (|2009L 1201 1[ ) accurately 
determine the UV-slope of LBGs at z > 3 as a function 
of UV-luminosity and redshift. For a star- forming stellar 
population, the UV-slope, parametrized by /3 {fx ex A^), 
is most sensitive to the dust content of the galaxy. In 
consequence, the authors argue that the dependence of P 
on UV-luminosity can be fully explained by an increase 
in the amount of dust that these galaxies contain as a 
they become brighter. If this is the case, a trend of red- 
der {Un — Vn)rest colors for UV-brightcr sources like the 
one observed here is also expected, since this color is 
also affected by dust extinction. The bottom panel of 
Figure [6] shows the {Un — Vn)rest color s corrected for 
dust e xtinction. These were derived as in lBouwens et al.l 
()2009f ). fro m the mean vs. UV-luminosity trend and 
the classic iMeurer et al.l () 19991 ) relation that relates the 
dust content with the UV-slope /3 according to: 

Ai6oo= 4.43 - 1.99/3 

whe re Aienn is t he ex tinction in magnitudes at 1600 A. 
A Calz etti et al.l ()2000l ) extinction curve is used to esti- 
mate the extinction at other wavelengths. 

As can be seen from the figure, the dust corrected {Un — 
Vn) color is essentially constant with UV-luminosity. 
This again suggests that the main origin for the correla- 
tion of P with the UV-luminosity is an increasing dust 
content at brighter magnitudes. Otherwise we would ex- 
pect some residual trend in the dust corrected color. This 
would be the case for example, if a strong age trend ex- 
isted with UV-luminosity. As will be shown later, a sim- 
ilar flattening of the corrected colors is observed at all 
redshifts in the — 21 < Mi 500 < — 18 range. 

Although both are flat, the trends of the corrected col- 
ors with UV-luminosity of the two samples are disjoint. 
The reasons for this are unclear but there are a couple of 
reasons that may help explain it. First, it is unclear that 
the mean trend used to derive dust corrections applies 
over such a large range of luminosities. This can help to 
explain that the corrected colors of the brightest sources 
are bluer than those of the faintest ones. However, it 



TABLE 2 

Photometry Summary of Stacked SED. 



Ref. mag Ners Nhudf B435 Veoe irrd zsbo Yos Yiob J125 Hieo [3.6] [4.5] 



Z775 z ~ 4 5— dropouts 

25 27 / 55 / 4 27.5 ± 7.2 224.0 ± 14.9 363.2 ± 18.8 384.8 ± 24.7 366.5 ± 22.4 ••• 405.5 ± 33.0 480.7 ± 53.4 784.5 ± 78.3 768.6 ± 117.2 

26 87 / 16 12 / 17 13.5 ± 2.4 103.1 ±3.8 150.1 ±4.1 149.1 ± 5.8 146.8 ± 5.9 130.6 ± 8.9 147.9 ± 5.3 161.8 ± 7.1 268.8 ± 17.1 280.0 ± 17.6 

27 107 / 2 25 / 34 8.4 ± 1.2 45.5 ± 1.7 58.9 ± 1.8 54.4 ± 2.4 55.0 ± 2.9 53.1 ±3.2 56.1 ±2.3 60.9 ± 3.4 83.8 ± 6.3 67.0 ± 7.5 

28 45 / 78 43 / 74 3.6 ± 1.0 24.4 ± 1.0 27.2 ± 1.0 26.0 ± 1.6 29.4 ± 3.3 22.4 ± 1.5 26.8 ± 1.5 29.0 ±2.1 41.4 ±4.9 43.7 ± 7.2 



Z850 z 5 y— dropouts 

25 4 / 9 / -6.2 ± 11.7 42.1 ±8.1 227.4 ± 40.8 267.8 ± 21.8 273.0 ± 22.6 ••• 281.9 ± 18.4 258.3 ± 19.3 681.9 ± 62.2 498.8 ± 59.9 

26 18 / 28 2 / 3 0.9 ± 3.0 15.9 ± 3.7 104.7 ± 11.9 125.6 ±8.1 114.5 ± 13.2 128.4 ± 17.4 126.3 ± 8.7 137.3 ± 11.1 337.9 ± 40.3 236.5 ± 41.5 

27 43 / 61 8 / 9 -2.9 ±2.1 8.3 ± 1.7 58.9 ±4.2 63.9 ± 3.0 70.0 ± 5.6 62.1 ±8.0 67.7 ±3.9 63.0 ± 4.6 116.3 ± 11.5 89.4 ± 10.6 

28 11 / 24 9 / 13 -0.8 ±1.9 3.5 ± 1.4 24.5 ± 3.0 26.8 ± 2.6 25.5 ± 6.2 19.0 ± 3.7 25.9 ± 3.9 26.5 ± 3.9 58.0 ± 15.4 33.8 ± 15.7 



Y z 6 i— dropouts 

26 3 / 6 3 / 3 1.2 ±3.8 -1.5 ±3.6 12.1 ±6.6 93.0 ± 14.7 102.2 ± 12.1 143.2 ± 23.8 116.4 ± 14.8 133.6 ± 26.3 278.7 ± 56.8 187.7 ±42.1 

27 11 / 17 2 / 3 -0.5 ±3.4 0.1 ± 2.4 14.2 ± 5.0 64.5 ± 6.5 60.2 ± 7.4 48.7 ±5.1 64.2 ± 5.9 68.4 ± 7.8 142.2 ± 27.0 111.9 ± 31.8 
28 4/4 6/13 -2.5 ±2.5 0.5 ±2.1 -0.7 ±3.6 30.3 ± 5.4 30.7 ±8.9 24.6 ± 3.4 33.0 ± 5.8 26.2 ± 6.3 61.3 ± 19.2 39.9 ± 20.0 



HiQo z ^ 7 2;— dropouts 

26 ••• ••• -3.2 ±2.3 2.8 ±2.3 -7.4 ±4.7 23.5 ± 5.4 63.7 ±7.7 ••• 129.7 ± 16.0 128.2 ± 17.0 262.0 ± 28.0 181.0 ± 37.0 

27 ••• ••• -2.4 ±0.8 -0.8 ±0.6 2.2 ± 1.8 14.6 ± 2.3 49.1 ±4.3 ••• 64.8 ±4.1 57.0 ± 3.4 107.0 ± 16.0 83.8 ± 25.0 

28 ••• ••• 1.2 ±0.8 0.4 ±0.5 -0.1 ±0.7 5.6 ±1.1 27.2 ± 1.8 ••• 29.6 ± 2.0 24.5 ± 1.6 45.1 ±9.5 39.3 ± 17.0 



Note. — The median SEDs of high— 2; galaxies in bins of observed UV- luminosity. The fluxes are in units of nJy. Nhudf and Ners are the numbers of sources with r eliable IRAQ 
photometry over the total available in the given bin. Only the ones with reliable photometry were considered in the stack. The 2; ~ 7 stacks were taken directly from lLabbe et all 
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TABLE 3 

Rest Frame Magnitude and Colors. 



(C^n-K).e.. = -0.026 xMi5oo-0.069 
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Fig. 6. — Top: Interpolated rest frame (Un — Vn) color as a 
function of M1500 for the 2; ~ 4 sources. Blue solid symbols and 
errorbars correspond to the colors derived for the new SEDs in this 
work. Similarly, we have derived interpolated {Un — Vn)rest colors 
for the SEDs presented in Lee et al. (2011), dark open symbols and 
errorbars. Although not very steep, a systematic trend to bluer 
(Un — Vn)rest colors with decreasing UV-luminosity is observed. 
The best fit trend of (Un - Vn)rest = -0.026 x M1500 - 0.069 is 
shown by the dashed line and the 95% confidence interval is shown 
by the gray area. The rest-frame (Un — Vn) colors were determined 
from best fits to the observed data based on the default template 
set from the code EAzY ( Brammer et al. 2008). We have used ideal 
narrow band Un and Vn filters that correspond to step functions 
of 100 A width, centered at rest-frame 3500 A and 5500 A respec- 
tively. The upper limit (2cr), arises from the limitations of the 
template set which only includes (Un — Vn) > O.S mags. The color 
of this SED is likely very close to this limit. Bottom: The same 
(Un — Vn)rest colors after being corrected for dust extinction. The 
corrections were derived from the median UV-slope (3 t rends from 
iBouwens et al. (2009, 2011) and the local IMeurer et al . (1999) re- 
lation between $ and dust extinction. The corrected colors of both 
samples show fiat trends with luminosity but seem to be disjoint. It 
is possible that the dust corrections are not equally adequate over 
the entire luminosity range. Another possible explanation is that 
the contamination rates are different between the two samples. Fi- 
nally, it should be noted that the UV-slope in the Lee et al. (2011) 
sample is only poorly sampled. Small offsets in the photometric 
calibration could bias the best fit models, from which the colors 
are estimated, to redder colors. 

does not explain the break between the two samples. Sec- 
ond there are differences between the selection of the two 
samples, in particular, it is not clear that the fraction of 
contaminants is comparable between the two samples. A 
larger fraction of contaminants in the brighter ground- 
based sample of lLee et al.l (|2Qlll ) as compared to the our 
data, could bias the brighter SEDs to smaller {Un — Vn) 
colors, resulting in excessive corrections for dust, which 
are derived from the mean trend on the UV with UV- 
luminosity. I^_should also be noted that in the case of the 
ILee et al.i ()2Qllf ) sample, the UV-slope is not very well 
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Note. — Rest frame M1500 estimated from the best fit power- 
law to the rest-frame UV photometry. Possibly contaminated 
bands were ignored in the fit (see open points in Figure [T]). 

sampled, with basically the / and J— bands driving the 
UV-slope (since the i^— band has been excluded). Even 
small systematic offsets of ~ 0.1 mag in the zero-point 
calibration between these two filters can bias the best 
fit models to redder slopes and in consequence, smaller 
UV-to-optical colors. 

5.2. The UV-to-optical Colors at Higher Redshift 

We now extend our determination of the SED of high-z 
galaxies to 2: ~ 5 and z ~ 6, and combine our determina- 
tions with th e stacked SEDs of so urces at 2: ~ 7 recently 
presented in iLabbe et al.l (|2QlQa[ ) . Figure [7] shows the 
SEDs at 4 < z < 7. The expected mean redshifts of these 
samples correspond to (z) = 3.8, (z) = 5.0, (z) = 5.9, and 
(z) = 6.9, with typical spread in redshift of Az±0.3. The 
samples at each redshift have been split in bins of 1 mag, 
using as reference the observed magnitude in the band 
closest to 1500 A (Table [2] contains the SEDs and infor- 
mation about the samples in each bin). Best fit BC03 
models with CSF and 0.2 Zq metallicity are overlaid for 
reference. As explained in the previous section, we have 
excluded from the fitting the bands that are potentially 
biased (marked by the open symbols). As can be seen in 
this Figure, all these SEDs are remarkably similar and 
show sizable UV-to-optical colors. 

The similarity of the SEDs can also be appreciated 
in Figure [HI where the SEDs from the different redshift 
samples have been grouped according to their approxi- 
mate M1500 luminosities. The SEDs in the three M1500 
bins have been re-normalized according to their M1500 lu- 
minosity and combined to produce over-sampled SEDs, 
which take advantage of the fact that the different filters 
probe different rest-frame wavelengths at the different 
redshifts. The combined SEDs show the typical flat UV 
slopes, but also rather fiat colors in the rest-frame op- 
tical, albeit with a relatively large scatter. This scatter 
can in part be a consequence of the normalization which 
is done in the far UV, which effectively acts as the pivot 



9 



1 1 — I I I I — T"^ — ^ — — I — p 1 — I I I I — I \ — ;; I ; I 




O O O O >^ >^ >^ <? ^ O O O O >^ <? ^ 



Observed Wavelength [//m] 



Fig. 7. — The stacked SEDs of galaxies in units of observed magnitudes (see also Table O. The x-axis shows wavelength and the filter 
that it corresponds to for reference. In the case of the optical bands, the errors were derived by bootstrap re-sampling the measured fiuxes. 
The individual uncertainties were set to a minimum of 5% to account for systematic uncertainties in the filter to filter absolute calibrations. 
The errors in the IRAC bands were derived by bootstrap re-sampling the individual images and repeating the pixel by pixel median-stack 
process. In all cases, these errors on the median should include both the image noise and the uncertainty coming from variations within the 
population included in the stack. Simple best fit BC03 models with CSF are also included for reference. In the fitting process the redshifts 
were fixed to the median redshift of each sample, i.e.: (z) = 3.8, (z) = 5.0, (z) = 5.8, and (z) = 6.9 respectively. Due to the intrinsic 
redshift distribution of each sample, the median fiuxes measured for bands near Lyman alpha or the Balmer break are likely biased and 
in consequence they were excluded from the fitting process (excluded bands are marked by the open symbols). The position of the most 
prominent possible emission lines in the rest-frame optical region of the SED are marked by the vertical dashed lines (assuming the mean 
redshift of each sample). At the different redshifts sampled here, different lines could be contaminating the continuum fiuxes measured by 
the Spitzer /IKAC filters. The overall shape and UV-to-optical colors of the SEDs, however, remain remarkably constant with redshift. 
Some contribution of nebular emission lines is possibly hinted by the small relative excess of [3.6] over [4.5] fiuxes observed at z > 5. 



point. Different emission lines that sometimes contribute 
flux in the IRAC bands, depending on the redshift, can 
also be contributing to this scatter. 

5.2.1. Emission lines in the rest- frame optical. 



We first consider the implications if the sizable red 
-^125 — [3.6] colors that we observe are assumed to be ex- 
clusively measuring the stellar continuum in these high 
redshift galaxies. The red color would arise from large 
Balmer breaks, and would be indicative of evolved stellar 
populations. As can be seen in Figure [9^1, the J125 — [3.6] 
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Fig. 8. — The SEDs of the combined samples in three bins of 
intrins ic U V-luminosity: M1500 ~ ~21, —20, and -19. As seen in 
Figure [To] the UV-to-optical colors of the SEDs depend weakly on 
redshift. Here we perform the exercise of combining all the photom- 
etry from SEDs at different redshifts into a single SED. At different 
redshifts, our set of HST/ACS+WFC3/IR+S'pifzer/IRAC filters 
probe different rest-frame wavelengths. The combined SEDs are 
plotted in terms of their magnitudes relative to M150O5 which is 
measured from the best power law fit to the UV continuum of each 
median SED. An offset magnitude of —21, —20, and —19 is added 
to reflect the approximate intrinsic magnitude of the SEDs that 
are being combined. 

colors characteristic of these sources would then require 
fairly large ages of ~ 1 Gyr if no dust is allowed in the 
models (the model tracks in the Figure are CSF models). 
Allowing some re ddening by dust (using in this case a 
ICalzetti et al.ll2QQQ law), reduces their ages considerably 
(Figure[9t). The amount of dust reddening is constrained 
however. Based on the UV conti nuum slopes typically 
observed in th ese galaxies ( Bouwen s et al.ll2QQ9f ) , and us- 
ing the z = Meurer et al. (1999) relation between this 
slope and the dust content, it has been shown that the 
reddening at z > 4 is typically Ay ^ 0.7 mag. The 
redde ning is likely to be even smaller at higher red- 
shifts (jBouwens et al.ll2QQ9L l2QTTh . With this added con- 
straint, the ages derived from SED fitting are generally 
estimated to be 300 - 400 Myr ( e.g., Yan et al. 2005; 
lEvles et a l. 2005': ' Stark et al.ll2009l : [Gonzalez et al.1120101 : 
iLabbe et a l. 2010a^. 

In our SEDs at 2: > 5 we also find that the ([3.6] — 
[4.5]) color is consistently blue, with typical values of 
~ —0.3 mags. This color is very hard to reproduce with 
a model with any star formation history that is based 
only on stellar continuum. The addition of dust only 
makes the comparison worse (Figure [9fc,d). It should be 
noted, however, that the uncertainties in this color are 
considerable and that in most cases the observed color 
is formally consistent with the stellar continuum only 
models. 

It has been claimed (Schaerer fc de Barrosl l2009ll2010f ) 
that at some of these redshifts, particularly z ^ 6 and 
7, the SEDs could be better fit by very young stellar 
population models (< 10 Myr old) with moderate dust 
content and very strong nebular emission lines. These 
lines would dominate the rest-frame optical fluxes. As 
the galaxies are forming stars, it is likely that they have 
nebular emission lines, and so we have assessed whether 
emission lines could be influencing our fluxes and fits. 
The observed position of the most prominent optical neb- 
ular emission lines at each redshift is marked in Figure [71 
by the dashed vertical lines (assuming the mean redshift 



of each sample). It can be seen that at the different red- 
shifts involved, some of the rest-frame optical emission 
li nes would li e withi n the IRAC [3.6] and [4.5] filters. 

lAtek et al.l (120111 ) showed recently that at z < 2.8 
there are sources that show emission lines with large 
enough equivalent widths to dominate the broadband 
fluxes in the optical. However, even for their sample 
of high equivalent width selected sources, such objects 
are not very common. For the more typical cases, the 
lines that they find in their sample would have moderate 
contributions to the IRAC broad bands, usually mak- 
ing them 0.2 — 0.3 mags brigh ter. At highe r redshifts, 
3.8 < ^ < 5, the work of Shim et al.l (j2011f ). based on 
spectroscopic redshifts and IRAC broadband photome- 
try, finds similar results. They find that the sources in 
this redshift range show an excess of flux in [3.6] over 
[4.5], which the authors think can be explained by Ha 
emission. The majority of the sources in their sample 
have ([3.6] — [4.5]) colors that are consistent with an Ha 
contribution of 0.2 — 0.3 mags to the [3.6] filter, indicating 
moderate contributions from emission lines. 

We consider the effects of a similar (moderate) con- 
tribution of emission lines to the colors of a 280 Myr 
old CSF model. For the purpose of this exercise we 
assume a Ha rest-frame equivalent width EWrest = 
300 A. At z ~ 6, this would increase the luminosity 
in the [4.5] filter by 0.23 mag. We only consider the 
strongest lines in the optical and assume the str eng th ra - 
tios presented by Anders & Fri tze-v. AlvenslebenI ()2003l ) 
for a 0.2 Z0 metallicity system. This corresponds to 
EWrest (OH, H/3, OIH) = (189, 105, 670) A, respectively. 
The open squares in Figure [9^1, b, correspond to the me- 
dian colors for a model like the one described (the median 
takes into account the expected redshift distribution of 
each sample). 

As can be seen in Figure [9^, the main effect of the 
contribution from the emission lines is that it allows for 
significantly younger ages to reproduce the J125 — [3.6] 
colors (by roughly a factor of two). Since this implies 
lower stellar continuum fluxes than previously assumed, 
this would also result in lower best fit stellar masses than 
previously estimated for galaxies at these redshifts (by a 
similar factor). Simultaneously, given the line ratios and 
the redshift distributions expected for these samples, a 
model like this can also help in reproducing the blue [3.6]- 
[4.5] colors of the SEDs at z > 5 (Figure Eb). While for- 
mally models with just stellar continua can fit the data, 
given the current large uncertainties, such colors are hard 
to obtain with stellar continuum only models. It is more 
likely that some degree of contribution of emission lines 
provides a better solution. 

The assessment presented here is necessarily very sim- 
plistic, given the current uncertainties, and is only meant 
to show the possible effects of emission lines on the colors 
of these galaxies. Better estimates of the contribution of 
these lines and its effects on the properties derived from 
SED fitting will depend on the actual strengths and ra- 
tios of the lines. Moderate improvement can be achieved 
if more precise redshifts are known for these sources, such 
data does not exist at this time for such faint high red- 
shift galaxies and so this remains an open question until 
improved spectroscopic capability becomes available. 
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Fig. 9. — The observed colors of the median SEDs. The measurements are shghtly offset in redshift for clarity, a) Theoretical tracks 
for a BC03 CSF model without dust are included. The different lines correspond to ages: 50, 140, 280, and 720 Myr. The 280 Myr 
old track is highlighted. For these dust-free models, the J125 — [3.6] color is more consistent with older galaxy populations in which 
there is a significant contribution of relatively evolved stars to the rest-frame optical light. The most likely age is lowered if flux from 
optical nebular emission lines is included in the appropriate IRAC band and redshift. As an example of the changes that can result 
from the inclusion of emission lines, the open squares (connected by the dashed line) show the effects of adding emission lines with 
EWrest(OII, H/3, OIII, ¥La) = (189, 105, 670, 300) A, respectively to the 280 Myr old track, b) Dust free tracks (for the same ages as in 9a) 
cannot reproduce the [3.6] — [4.5] ~ —0.3 colors observed at 2; > 5. The inclusion of optical nebular emission lines, shown by the open 
squares connected by the dashed line can help reproduce such blue colors, c) As in a) but only the 280 Myr old track is shown with 
varying amounts of dust added. The tracks correspond, from bottom to top, to Ay = 0.0, 0.2, 0.4, and 0.6 mag using a Calzetti et al. (2000) 
extinction law. As expected, redder colors can be obtained for younger models if dust is allowed. Our independent analysis of the UV 
continuum slope (Bouwens et al. 2011), constrains the amount of dust to Ay < 0.7 mag for z > 4. With this added constraint, the ages 
of high-z galaxies are generally estimated to be > 400 Myr when the emission line contribution is taken to be zero, d) As is clear here, the 
addition of dust makes it even harder for the stellar continuum only models to reproduce the blue [3.6] — [4.5] colors observed in the SEDs. 
This suggests that an optical nebular emission line component to the total flux may be needed. 



5.2.2. UV-to-optical color vs. luminosity. 

Regarding the trend of redder UV-to-optical colors for 
brighter UV-luminosities presented in Figure [6] for the 
z ~ 4 sample, it can also be observed for the other sam- 
ples. In fact, within the uncertainties in the color deter- 
mination, the J125 — [3.6] vs Ml 500 trend is roughly the 
same at z ~ 5 and 6 (Figure [TQ| top). This is particularly 
remarkable in view of the fact that this color is an ob- 



served color, not a rest-frame color, and in consequence 
is probing different wavelengths for the different samples. 
This shows again how similar and flat these SEDs are. 
If all the SEDs are considered simultaneously, a trend 
of J125 - [3.6] = -0.12(±0.07) X Mi5oo,AB + const, is 
found. For simple CSF models, this trend implies a vari- 
ation in stellar M/L ratio. Gonzalez et al. (2011) find 
that at z ~ 4, this variation is of a factor ~ x5 between 
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Fig. 10. — Top: UV-to-optical color as a function of Mi 500 • Blue 
symbols correspond to the J — [3.6] colors of the z ~ 4 SEDs; green: 
2; ~ 5; red: z ~ 6; and black: z ^ 7. Despite the sizable scatter, a 
systematic trend to redder UV-to-optical colors with increasing lu- 
minosity is observed at all redshifts. At the same time, there is no 
clear indication of a systematic trend with redshift. The dashed 
line is a best fit to all the points (all redshifts simultaneously): 
-^125 — [3.6] = — 0.11(±0.07) X Mi5oo,AS + const, and the gray 
area corresponds to the 95% confidence interval. This is different 
to Figure [6] in that this Figure shows the directly observed colors 
rather than interpolated rest-frame colors. In fact it is remarkable 
that this observed color, which samples different wavelength re- 
gions of the SED at the different redshift, remains fairly constant. 
This is probably an indication that the rest-frame optical colors 
measured by IRAC are fairly fiat. Center: The interpolated rest 
frame {Un — Vn) colors (Section 5.1). The best fit trend to all the 
data is shown by the brown dashed line. For comparison, the black 
dashed line shows the best fit to the z ~ 4 data only, as shown in 
Figure [6] this trend is slightly fiatter. Bottom: The (Un — Vn)rest 
colors after being corrected by dust extinction followi ng a simple 
prescription based on the UV-slopes and the local M eurer et al.l 
|"1999) relation. Except for the z ^ 7 SEDs (which have a much 
steeper (3 vs. UV-luminosity relation), the dust corrected colors of 
all SEDs seem to show a fiat trend (although the scatter is large). 
This indicates that a change in dust is enough to explain the change 
in colors, both in the UV and the UV-to-optical. In particular, no 
age dependence on luminosity seems to be required. 



^1500 = — 18 and M1500 = —21. So, despite the steep 
UV-luminosity functions characteristic at z > 4, the con- 
tribution of the faintest sources to the total stellar mass 
density is more modest than their contribution to the 
star formation rate density. 

Similar to the analysis performed at z ~ 4, we have 
derived rest-frame iUn — Vn) colors for the median SEDs 
at all redshifts. These colors are shown in Figure [TOl 
The trend with UV-luminosity has large scatter and low 
significance but is systematic in the sense that brighter 
sources are redder in {Un — Vn)rest' A fit to the rest-frame 
colors as a function of M1500 of all redshifts simultane- 
ously, results in (Un — Vn) = —0.04 x M1500 — 0.37, shown 
as the brown dashed line in Figure [10] {center panel). 
This is slightly steeper than the —0.03 slope found for 
the z ~ 4 data only, also plotted for reference as the 
black dashed line in the same Figure. The latt er trend 
included the bright sample of Lee et al.l (12011') and is 
consistent within the uncertainties with the fit to all the 
points. 

Usin g the median trends determined bv lBouwens et all 
(|2011f ) regarding the j3 slopes vs. UV-luminosity for 4 < 
z < 1 LBGs, we corrected th e (U — V)rest c olors for 
dust extinction using the Meu rer et al.l (j 19991 ) relation 
(Figure [TOt bottom). These corrected colors show large 
scatter around {U — V) rest, corrected ~ 0.4 but there is no 
indication of a residual trend with UV-luminosity except 
for the z 7 sample (black symbols). This may indicate 
again that at 4 < z < 6, the evolution can be fully 
explained by a change in dust content and there is no 
strong evolution in the ages of these sources as a function 
of UV-luminosity. This is consistent with the report by 
(jStark et al.ll2009l ) based on stellar population modeling 
of slightly brighter LBGs in this redshift range. It is also 
consistent with the prediction s from the hydrodynamical 
simulations of lFinlator et al.l ()2011f ). which find no trend 
of age with UV-luminosity. 

A.i z ^ 1 there seems to be a residual trend in 
(Un - Vn)rest,corrected ^ 0.4 VS. UV-lumiuosity in the 
sense that brighter sources have systematically bluer in- 
trinsic (Un — Vn) colors. In principle, this could indicate 
that the brighter sources are actually younger than the 
fainter sources, which would be difficult to imagine. Al- 
ternatively, it is possible that the simple dust corrections 
(calibrated at 2: = 0) are not adequate at these high red- 
shifts, due to, for example, changes in the IMF, metal- 
licity, or escape fraction, all of which could alter the UV- 
slope f3 for a given dust content. It should also be noted 
that the /3 vs. UV-luminosity at z ~ 7 is considerably 
steeper than at the other redshifts and so the corrections 
to the {Un — Vn)rest VS. UV-lumiuosity trend are more 
extreme. Nevertheless, the uncertainties are sufficiently 
large that also the z ^ 7 trend is fully consistent with 
being constant. 

6. SUMMARY 

Determining the rest-frame UV-to-optical colors of in- 
dividual UV-faint galaxies at z > 4 is challenging. At 
such redshifts the rest-frame optical lies at A > 2.5 mi- 
crons in the mid-IR and so it is hard to measure at the 
extremely faint magnitudes typical of galaxies in the ffrst 
1.5 Gyr. To gain insight into the typical spectral proper- 
ties of z ~ 4 — 6 galaxies, we have taken advantage of the 
large samples of 2; ~ 4 — 6 sources found in the deep HST 
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ACS and WFC3/IR images of the ERS and HUDF fields. 
The HST data give us rest-frame UV fluxes. We then use 
Spitzer /IUAC [3.6] and [4.5] micron data from the deep 
GOODS survey of these fields to get rest-frame optical 
fluxes by stacking the IRAC images and determining me- 
dian flux values. This allows us to determine the SEDs 
of these galaxies covering both the rest-frame UV and 
the optical. These SEDs represent a flrst comprehensive 
study of the rest-frame UV-to-optical properties of 2: > 4 
galaxies as a function of both redshift and UV-luminosity 
down to very faint limits. 

At z > 4 we have also combin ed our fai n t SED s with 
the L > L* stacks presented in iLee et "all ()2Qlll ). This 
allows us to assess the colors of 2; > 4 sources over an 
unprecedentedly large range of luminosities. Our main 
flndings are as follow: 

• At z > 4 the {Un — Vn) rest-frame color (inter- 
polated from the data using a fitted SED from 
a set of templates is of {Un — Vn) ~ 0.5 mags. 
There is a shallow but systematic trend of red- 
der colors for brighter UV-luminosities over the 
range —23 < Mi5oo,ab ^ —17.5 (Figure [6l section 
5.1). A linear fit to this trend is {Un — Vn)rest = 
-0.03xMi5oo -0.07 

• The SEDs of star forming galaxies at 4 < 2: < 7 are 
remarkably similar at all luminosities from — 23 < 
Misoo,AB ^ —17.5, showing fairly fiat rest-frame 
UV and rest-frame optical colors (section 5.2 and 
Figures \7\ [51 and [TQ|) . At a given redshift, there are 
weak indications of a subtle trend for redder colors 
for brighter sources. A simple fit to the data at all 
redshifts simultaneously results in {Un — Vn)rest = 
-0.04xMi5oo-0.37 (FigureHO] center. Section 4.2). 

• The UV-to-optical color, as measured by the ob- 
served J125 — [3.6], color remains fairly constant 
with redshift, despite the Spitzer /IRAC bands 
probing different wavelength regions of the SED 
(Figure [To] ^op). This suggests that the optical col- 
ors are fairly fiat at all redshifts z ~ 4 — 7, although 
there is significant scatter, which can be caused by 
optical emission line contamination to the IRAC 
filters. 

• The z > 5 SEDs show consistently blue [3.6] — 
[4.5] ~ —0.3 colors. This is hard to reproduce with 
models that only include stellar continuum, al- 
though it is still formally possible given the current 
uncertainties. Nonetheless, it is probably better ex- 



plained with moderate fiux contributions from op- 
tical nebular emission lines in the two IRAC bands. 
Including some contribution from emission lines 
leads to lower best fit ages and lower stellar masses 
than previously estimated for galaxies at these red- 
shifts. Changes of a factor 2 could well result, but 
we would caution that detailed assessments are not 
yet possible with the current data. 

• We derive dust-corrected {U — V) colors using 
a simple dust correction based on the median 

UV-slope p v s. U V-lum inosity relation from 

[Bouwens et al.[ (j2009[ . [2011[ ) and the [Meurer et all 
([T999 ) relation (Figure ED\ bottom). The dust- 
corrected UV-to-optical colors at 4 < 2: < 6 are 
{U — V) ~ 0.4 mags. They do not appear to depend 
on UV-luminosity. This may suggest that moder- 
ate changes in the dust content of galaxies can ex- 
plain the color dependency observed both in the 
UV-slopes and the UV-to-optical colors. In partic- 
ular, it seems that no age evolution is required to 
match this color dependency. 

The recent deep and ultra-deep HST ACS and 
WFC3/IR imaging programs over the GOODS-S field 
have resulted in large samples of z ~ 4 galaxies. Tak- 
ing advantage of these large samples, along with a care- 
ful stacking analysis, we present a first comprehensive 
study of the typical SEDs over the rest-frame UV and 
the rest-frame optical for faint star forming galaxies in 
the 4 < z < 7 redshift range. The stacked SEDs allow 
us to study the UV-to-optical colors of high-z star form- 
ing galaxies down to very low luminosities. We find a 
weak trend to bluer color at fainter luminosities. Inter- 
estingly, these stacked SEDs also show a remarkable sim- 
ilarity across redshift in a period of considerable stellar 
mass growth from about 0.8 Gyr to 1.5 Gyr, suggesting 
a smooth self-similar mode of evolution at 2: ~ 7 — 4. 
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APPENDIX 
A. DEPTH OF THE STACKS 

To quantify the rest-frame optical properties of very faint high-redshift galaxies, we stack the IRAC images for a 
large number of faint galaxies, after removing the flux contribution from nearby neighbors (section 4). While we do not 
expect large systematics in the photometry of individual sources from our deblending procedure to remove neighboring 
sources, it is possible that modest systematics could arise when stacking a large number of sources. The role of this 
appendix is to determine how large such systematics might be. Fortunately we are able to demonstrate that the fluxes 
that we measure in the stacked IRAC images are unbiased and that their significance is properly determined. 

We started by selecting 200 empty areas in the HUDF field. We select them based on the segmentation maps derived 
from the high resolution HST images. We processed each of the selected empty areas as if they were the position of 
one of the real sources in our catalog, i.e. we ran them through the de-blending code to remove the flux from nearby 
un-associated sources. To determine the rms for a given Ngtacked number of stacked stamps, we randomly draw (with 
replacement) Ngtacked stamps from the 200 empty areas, median combine them and perform aperture photometry in 
the median stack in the same fashion as for our real stacks. We repeat the drawing 300 times and meas ure the r ms 
of these 300 trials. As can be seen in the upper panel of Figure [TT| the limiting flux decreases as ex l/V^stacked^ as 
expected for background-limited noise. In the case of the stacks that contain real sources, the rms is larger than the 
image noise because it is also affected by the intrinsic variation in the distribution of fluxes that are being stacked. 

As can be seen in the lower panel, the flux of the median stacks is very close to the expected zero flux. 
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Fig. 11. — Upper: The la limiting flux (in /iJy) as a function of the number of IRAC stamps stacked. AH the stamps have the same 
integration time of ~ 23.4 h. The estimates (open symbols) come from stacking large numbers of empty areas of the GOODS images. Each 
empty area was cleaned of neighboring fluxes using the code described in the text. The noise decreases oc with n being the number of 
empty regions stacked. The fllled symbols are the actual RMS determined for each of our stacked stamps: blue for the 2; ~ 4 stamps, green 
2; ~ 5, and red z ~ 6 (see Figures [21 [S] and[4|). These RMS values are > than the ones derived for empty regions. This is expected because 
the uncertainties in the stack include both the image noise and intrinsic variations in the population being stacked. Lower: The flux of the 
median combined stacks of cleaned empty regions as a function of the number of stamps combined. As expected, this is very close to zero. 



